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We discuss massive quark effects in the endpoint region x — > 1 of inclusive deep inelastic scattering 
(DIS), where the hadronic final state is collimated and thus represents a jet. In this regime heavy 
quark pairs are generated via secondary radiation, i.e. due to a gluon splitting in light quark initiated 
contributions starting at 0(a 2 ) in the fixed-order expansion. Based on the factorization framework 
for DIS in the endpoint region for massless quarks in Soft Collinear Effective Theory (SCET), 
which we also scrutinize in this work, we construct a variable flavor number scheme that deals 
with arbitrary hierarchies between the mass scale and the kinematic scales. The scheme exhibits a 
continuous behavior between the massless limit for very light quarks and the decoupling limit for very 
heavy quarks. It entails threshold matching corrections, arising from all gauge invariant factorization 
components at the mass scale, which are related to each other via consistency conditions. This is 
explicitly demonstrated by recalculating the known threshold correction for the parton distribution 
function at 0(q. 2 CfTf) within SCET. The latter contains large rapidity logarithms ~ ln(l — x) 
that can be summed by exponentiation. Their coefficients are universal, which can be used to 
obtain potentially relevant higher order results for generic threshold corrections at colliders from 
computations in DIS. In particular, we extract the 0(a 3 ) threshold correction multiplied by a single 
rapidity logarithm from results obtained earlier. 


I. INTRODUCTION 

For multiple collider processes precision phenomenol¬ 
ogy requires a thorough understanding of quark mass ef¬ 
fects to meet future experimental data with high statis¬ 
tics and low systematic uncertainties. One of the chal¬ 
lenges hereby is that the quark mass represents a scale 
that can be parametrically much different than the typ¬ 
ical hard scattering scale, which can give rise to large 
logarithms in fixed order perturbation theory. At the 
Large Hadron Collider (LHC) this concerns e.g. bottom 
mass effects in Higgs production via gluon fusion or in 
association with bottom quark jets. 

One example where the treatment of massive quark 
effects raised a lot of interest is deep-inelastic scattering 
(DIS), the benchmark process for the extraction of parton 
distribution functions (PDFs). These are one main input 
for the analysis of all processes at hadron colliders, such 
that precise predictions including the effects of the charm 
and bottom quark masses are necessary. A first system¬ 
atic approach to incorporate heavy quarks with arbitrary 
masses with respect to the other relevant scales has been 
provided in Refs. 0,1 , which laid the basis of a variable 
flavor number scheme (VFNS) for inclusive processes in 
hadron collisions. The method is founded on the separa¬ 
tion of close-to-mass-shell modes and offshell fluctuations 
and is thus in the spirit of effective field theory (EFT) fac¬ 
torization, see e.g. Ref. Q. Nowadays, different schemes 
have been developed to cope with this challenge, which 
are mainly based on this approach, but they differ in their 
detailed implementation concerning formally subleading 
contributions (see Ref. jdj for a short overview). 

In this work we discuss the endpoint region of DIS, i.e. 
x 1, where the final state becomes a single jet with 


invariant mass — % -C Q, where q 2 = —Q 2 denotes 
the hard momentum transfer. Although having a limited 
phenomenological impact, the regime 1 — x <C 1 provides 
a simple and instructive example of how to incorporate 
quark mass effects in differential distributions at hadron 
colliders with multiple kinematic scales. We construct 
a VFNS, which is in the same spirit as the well-known 
VFNS in the classical operator product expansion (OPE) 
region 1 — x ~ 0(1) 0[|j and exhibits similar main char¬ 
acteristics. These include (i) the resununation of all large 
logarithms, also those involving the quark mass, (ii) the 
correct limiting behavior of the perturbative structures, 
i.e. the hard matching coefficient at the scale Q and the 
jet function at the scale Qy/1 — x, for very small and large 
masses and (iii) a continuous description for arbitrary hi¬ 
erarchies of the dynamic scales with respect to the mass 
keeping the full mass dependence of the singular terms 
(i.e. at leading order in the power counting). The lat¬ 
ter is in particular relevant, since for a single value of the 
hard momentum transfer Q or the Bjorken variable x the 
hierarchies can change significantly when scanning over 
the respective other variable. 

We will see that for x —> 1 quark mass effects arise 
(mainly) via light quark initiated contributions, where 
the heavy quarks are produced via secondary radiation 
through the splitting of an additionally emitted virtual 
gluon starting at 0(a 2 ). Our treatment of these sec¬ 
ondary massive quark corrections relies on the setup de¬ 
veloped in Refs. i5j,@]. In these papers a VFNS for event 
shape distributions, specifically for thrust, in the dijet re¬ 
gion for e"*"e _ -collisions was constructed in the framework 
of SCET. It was shown that in a strict EFT interpreta¬ 
tion one has to introduce additional degrees of freedom to 
the existing massless SCET modes, namely collinear and 
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soft mass modes. These adopt the scaling of the corre¬ 
sponding massless modes if the mass is below the typical 
invariant mass scale of the massless collinear or ultrasoft 
modes, respectively, but contain in addition fluctuations 
around the mass-shell which have to be considered when 
the massive quark is integrated out. The mass mode 
picture leads to the emergence of different EFTs, which 
implies the strict guideline of having the massive quark 
modes either as fluctuating fields or excluded completed 
(i.e. integrated out). 

In Ref. @ also an alternative formulation of the VFNS 
is described which does not rely on different EFTs, 
but is only based on the massless factorization theorem 
with different renormalization conditions for the massive 
quark corrections to the matrix elements according to 
the hierarchy between the mass and the other involved 
scales. This renormalization procedure is similar to the 
one for the strong coupling in the presence of massive 
quarks 0 and is also the underlying idea for the for¬ 
mulation of the VFNS for DIS in the OPE region 0. 
The interpretation in terms of different renormalization 
conditions is convenient since it automatically takes into 
account any possible power corrections between the hard 
or jet scale and the mass scale arising in the EFT picture 
when the hierarchy between the mass scale and the re¬ 
spective kinematic scale is marginal!]] Thus it provides a 
continuous description of the cross section for arbitrary 
masses by construction, whereas in a strict EFT picture 
(which would enforce expansions) the transitions between 
different hierarchical scenarios would have to be adapted 
by nonsingular corrections, which concerns in particular 
the real radiation thresholds. In the following we will 
therefore discuss the VFNS only in the formulation rely¬ 
ing on renormalization conditions. 

The renormalization conditions with respect to the 
massive quark corrections which we are going to impose 
are either the MS prescription or an on-shell (OS) pre¬ 
scription!]] The common use of the MS prescription has 
the feature that the ni massless quarks and the massive 
flavor both contribute to the renormalization group (RG) 
evolution in the same way corresponding to an (n; + 1) 
running flavor scheme. The OS prescription is defined 
by the condition that the massive quark corrections van¬ 
ish for invariant mass scales much smaller than the quark 
mass and also subtracts finite and scale-dependent contri¬ 
butions such that the massive flavor does not lead to any 
contribution in the RG evolution, implying only rq run¬ 
ning flavors. This concerns the matrix elements, i.e. the 
current, the jet function and the PDFs, as well as the 
strong coupling a s . The MS prescription is appropri¬ 
ate to cover the situation where the quark mass becomes 


1 For all possible scale hierarchies covered within our VFNS ap¬ 
proach a strict EFT factorization can be constructed that agrees 
with our result up to nonsingular corrections. 

2 For the purely massless quark corrections we always use the MS 

scheme. 


small (where appropriate means that no large mass loga¬ 
rithms arise in this limit) and leads to expressions which 
give the known results for massless quarks in the limit 
m —> 0. The OS prescription is suitable to cover the 
decoupling limit, such that the effects of the massive 
quark vanish in the infinite mass limit. The differences 
of the renormalized quantities with respect to both of 
these renormalization prescriptions constitute matching 
factors, also called threshold corrections. Since the hard 
matching coefficient, the jet function and the PDFs are 
independent and in principle not exclusively tied to any 
particular factorization theorem, these factors represent 
also universal ingredients that appear in a similar way 
for the description of different processes. Here we will 
emphasize universal features of the threshold corrections 
and establish the connections to some of the results an¬ 
ticipated in Ref. Q. 

The outline of this paper is as follows: In Sec. lU we 
set up the notation and display the massless factorization 
theorem for the structure functions in DIS for 1 — x -C 1. 
Here we do not require to be in the kinematic region 
with the scaling 1 — x ~ Aqcd /Q as frequently adopted 
in the literature. We show that the factorization the¬ 
orem has the same form in the complete endpoint re¬ 
gion 1 — x > Aqqd /Q using the proper mode setup for 
1 — x Aqcd /Q- We also explain that massive quark 
effects can only arise via secondary radiation, and we 
show in Sec. m how to incorporate them consistently by 
setting up a VFNS for any gauge invariant component 
of the factorization theorem. For definiteness we discuss 
in Sec. m practical implementations of the VFNS for 
various hierarchies between the kinematic scales and the 
mass scale. Here we also consider different choices of the 
final renormalization scale in the factorization theorem 
that lead to consistency conditions between the thresh¬ 
old correction factors involved in the RG running of the 
corresponding matrix elements. In Sec. 0 we explicitly 
calculate the PDF threshold correction at 0(a 1 2 ) in the 
large x limit in the effective theory and show that our re¬ 
sult is consistent both with the one obtained in classical 
DIS 0 expanded for x — > 1 and with the result for the 
jet and hard function threshold corrections performed at 
0(a 2 ) in Ref. 0. In the endpoint region the threshold 
corrections contain large logarithms ~ a 2 s log related to 
the separation of the collinear and soft mass modes in 
rapidity, whose resummation we carry out explicitly via 
the rapidity RGE Haim. Based on the considerations 
in Secs. II V I and M we display also the explicit expressions 
for the threshold corrections up to the required order for 
a full N 3 LL analysis, which includes rapidity logarithms 
~ a 2 log, ~ a 3 log and ~ af log 2 . These results represent 
universal ingredients useful for various collider processes 
that involve PDFs in the endpoint region, jet functions 
and a hard function related to the one appearing in DIS. 
Finally, in Sec. IVIl we conclude. For the sake of compari¬ 
son, we provide in Appendix O the results for secondary 
massive quark corrections in the OPE region known from 
Ref. 0. In Appendix [Bl we show that for x —> 1 our re- 
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suits are in agreement with them. 


II. MASSLESS FACTORIZATION THEOREM 
FOR DIS IN THE ENDPOINT REGION 


with the kinematic constraint 0 < x < 1. We will work in 
the Breit frame, where q^ does not have an energy com¬ 
ponent and the initial state proton is n-collinear. Ne¬ 
glecting the proton mass the relevant momenta in the 
Breit frame in terms of lightcone coordinates read 


Before discussing quark mass effects we briefly describe 
the kinematic setup and the factorization theorem for 
DIS in the endpoint region 1 — x -C 1. Here we display 
the mode setup, highlight the relevant steps for its deriva¬ 
tion specifically for the hierarchy 1 — x Aqcd /Q and 
show that it can be readily combined with the commonly 
considered scaling 1 — x ~ Aqcd/Q- 


A. Kinematics of DIS 


_ Q n M . 
H 2 2 

P£ = 
x 2 


p “ = 


<9(1 - x) 

--- -nr 

2x 


( 2 ) 


In the endpoint region the hadronic final state is an n- 
collinear jet with an invariant mass P\ ~ Q 2 (l — x) -C 

Q 2 . 

The differential cross section for DIS can be decom¬ 
posed in terms of a leptonic and a hadronic tensor. The 
latter is defined by 


In the following we consider the scattering of an elec¬ 
tron off a proton via photon exchange. We denote the 
proton momentum by P M , the momentum of the incom¬ 
ing (outgoing) electron by k^ ( k '^), the incoming mo¬ 
mentum of the virtual photon by q^ = k!^ — k^ with 
spacelike invariant mass q 2 = —Q 2 <0 and the momen¬ 
tum of the outgoing hadronic final state X by The 
Lorentz invariant Bjorken scaling variable x is defined by 


i J d 3 4 ze iqz {P\T[J^(z)J v (0)]\P) , 

(3) 

with | P) denoting the initial proton state and the current 
e 2 . qi^^qi{z) summed over all quark flavors 
qi with corresponding electric charges e qi . We will just 
deal with unpolarized DIS, so that a spin average is al¬ 
ways implied. Using current conservation, which implies 
q^W^ = 0, one can decompose the hadronic tensor for 
the parity conserving vector current into the two struc¬ 
ture functions Fi(x,Q 2 ) and F 2 (x, Q 2 ), 


W ^( p ,q) = —Im 


W^(P,q) 



F 1 (x,Q) 


- g ^F 1 (x,Q) + - 



withg^ = g^— l/2(n M n 1/ +n M n 1 '). Here the longitudinal 
structure function Fl(x,Q) reads 

F l (x, Q) = F 2 (x, Q) - 2xF 1 (x, Q ), (5) 

in terms of Fi(x, Q) and F 2 {x, Q). These structure func¬ 
tions contain physics at different invariant mass scales 
and thus must be factorized to resum the corresponding 
large logarithms. 


B. Factorization Setup 

In this section we briefly discuss the derivation of 
the factorization theorem for inclusive DIS for massless 
quarks in the endpoint region 1 — x <C 1 and set up 
the notation employed for the rest of the paper. The 
factorization can be performed in a multi-step matching 


+ 


1 

~P~q 





n v \ 

Y jF L (x,Q). 


F 2 (x,Q) 


(4) 


procedure and has been carried out already a number 
of times HM3- However, here we focus on the proper 
mode setup in the limit Q( 1 — x) Aqcd and explain 
why the factorization theorem adopts the same form as 
for Q(1 — x) ~ Aqcd- Although this fact has been al¬ 
ready stated in several papers (e.g. Refs. [ami) , we 
believe that it is worthwhile to give a short derivation 
using our mode setup@ 

The relevant modes are displayed in Fig. [T] The n- 
collinear modes describing the initial state proton in 
the Breit frame always have the same scaling = 
(n-pn,fi-pn,Pn) ~ (Q, A 2 qcb /Q, Aqcd)- The final state 


3 We disagree with the mode setup in Ref. fl3tl which assumes 

nonperturbative messenger modes for the beam remnants at the 
invariant mass scale Aqcd\/1 — x <C Aqcd> while Ref. [Hjl never 
explicitly displays the scaling of the modes. 
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FIG. 1. Relevant momentum modes for inclusive DIS in the 
endpoint region x —> 1 with 1 — x Aqcd/Q- 


is strongly collimated for x —> 1 with a large momentum 
Q and an invariant mass Q\J 1 — x and is thus described 
by n-collinear modes scaling as p£ ~ Q( 1 — x, 1, y/1 — x). 
The kinematics in the Breit frame prohibits the appear¬ 
ance of a final n-collinear state, as can be seen from 
Eq. This has the important consequence that the n- 
collincar sector just enters the factorization theorem via 
a component which is local both in label space as well as 
in the residual coordinate, as has been also pointed out 
in Ref. 0. The remaining relevant low-energy modes 
contribute to the measurement of x or equivalently to 
the squared invariant mass ~ Q 2 ( 1 — a;) via a compo¬ 
nent n ■ p ~ Q(1 — x) (i.e. they have to lie on the ver¬ 
tical line below the n-coll. modes in Fig. [TJ. In fact 
all such modes give vanishing contributions in pertur¬ 
bation theory, since no physical scale is associated with 
the other momentum components which results in scale¬ 
less integrals. This holds in particular also for ultra- 
soft modes scalin g as Q (1 — x, 1 — x, 1 — x) as stated 
e.g. in Refs, [l2j, lH, H3- Thus any additional rele¬ 
vant modes can only be nonperturbative and scale like 
Vcs ~ (Q(l-z), AqcdAK 1- ®). a QCd)- They encode in¬ 
terference effects between soft initial and final state radia¬ 
tion. Note that in contrast to the case 1 — x ~ Aqcd/Q, 
where the corresponding modes adopt the soft scaling 
p £ ~ Aqcd( 1, 1,1), these modes are now also boosted in 
the Breit frame, and therefore referred to as collinear-soft 
(csoft) modes. They are separated by the rapidity factor 
(1 — x) from the h-collinear modes. These types of modes 
have recently received some attention in the context of 
multidifferential cross sections and have been incorpo¬ 


rated systematically into a modified version of SCET, 
called SCET + pjl[20j. We will discuss here the DIS fac¬ 
torization theorem in the same spirit using a multistage 
matching procedure. However, our case is simpler, since 
no relevant softer mode is present with which the csoft 
mode can potentially interact 0 

To derive the factorization theorem we employ the mul¬ 
tistep matching procedure sketched in Fig. [2] We empha¬ 
size that other ways to order some of the matching steps 
are possible and one may even work in a framework con¬ 
taining the csoft modes from the beginning. Here we first 
match the QCD current to the usual SCET I current 

^QCD * ^SCETI = C{Q, p) XnYn'y^YjiXn , (6) 


with Xn = wkn and Xn = where W n and W n 

denote collinear Wilson lines. The ultrasoft Wilson lines 
Y n and Ik appear in the current after the BPS field redef¬ 
inition [21] that disentangles the collinear and ultrasoft 


sectors in the Lagrangian. SCET I describes collinear 
fluctuations at the invariant mass scale Q\J 1 — x with 
residual momenta of order Q(1 — x). When lowering the 
virtualities the n-collinear final state modes need to be 
integrated out, and the ultrasoft sector is being resolved 
with momentum components ~ Q(l—x) becoming labels. 
We call the corresponding theory with csoft modes of 
virtuality > Aqcd and h-collinear modes of virtuality > 
Aqcd/ y/1 — x Aqcd SCET 1+0 Here Aqcd/\/1 — x 
is the invariant mass scale of the collinear modes at 
which they can still interact with the csoft modes via 
the momentum component h • pn ~ Aq CD /(Q(1 - a;)). 
The matching coefficient between the two theories is the 
quark jet function a vacuum correlator of the 

hard collinear fields in SCET I describing the production 
rate of an inclusive jet with invariant mass s. It is defined 
in terms of the n-collinear fields as 


J(Q r niH) = 


-1 


2 ttN c Q 


x Im 


d 4 


ze-- z ( 0|T{* n , Q (0)^Xn(*U|0> 


(7) 


where the invariant mass is rjj ~ Qr+ and Xn,Q = 
Sn-r,Q Xn with V M denoting the label momentum opera¬ 
tor. All color and spin indices are traced implicitly. Here 
and in the following all expressions are only given for ini¬ 
tial state quarks. For antiquarks the corresponding ex¬ 
pressions are completely analogous. The matrix element 
in SCET I + is the quark PDF given by 



■ p)V^nXn, Q \P) ■ 

( 8 ) 


4 Here only the separation in rapidity matters, i.e. the ’’softness” 
of the csoft mode with respect to the n-collinear mode, while the 
boost of the csoft mode in the Breit frame is actually irrelevant. 

5 We do not consider ultrasoft modes in SCET I_|_ since they do 
not contribute to any measurement as stated before. 
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FIG. 2. Schematic picture of the multistage matching procedure for 1 — x Aqcd /Q described in the text. 


Here V„ and X n are Wilson lines of the label and small 
component of the csoft fields required by gauge invari¬ 
ance, written in the notation of Refs. [19j, 12010 In our 
matching procedure they originate directly from boosting 
the ultrasoft Wilson lines in SCET I 


be written a^] 

fiq/pi 1 -z,n) = Q j (Mg q/P (Q(1 - z) - £, n) S c (£, p ), 

( 11 ) 


yt y -1 = p 

1 n 7 v n x 


exp 


Y n X n = P exp 


pOO 

ig / dsn ■ A cs (sn 11 + x^) 
Jo 

ig / d s fi ■ A cs (sn^ + x 1 *) 

J — OO 


( 9 ) 

( 10 ) 


The change of the modes in the soft sector involves in 
principle an additional matching coefficient. However, 
since the ultrasoft sector in SCET I and the csoft sector 
in SCET 1+ are fully decoupled from any other sector 
and this matching only involves a single sector in both 
theories, for which we can use the QCD Lagrangian, the 
form of Eqs. © and m makes clear that no nontriv¬ 
ial matching coefficient is generated. As a last step the 
fluctuations at the invariant mass scale Aqcd/V1 — x 
need to be integrated out to describe the physics at the 
scale ~ Aqcd, the invariant mass of the initial state pro¬ 
ton, in terms of the final EFT, which we call SCET II + . 
Since the interactions between the n-collinear and csoft 
sectors are already fully disentangled there is no match¬ 
ing coefficient originating from the collinear sector and 
the transition from SCET 1+ to SCET II + , where these 
two types of modes cannot interact with each other, and 
can be achieved simply by lowering the virtuality of the 
n-collinear modes, in analogy to the matching between 
SCET I and SCET II [H]. The PDF in SCET 11+ can 


where g q /p(i , g) denotes a local collinear matrix element 

g q /p{l,V) = ^{P\Xn(0)^Xn, Q mP)m, ( 12 ) 

and S c (£. ji) denotes a vacuum expectation value in terms 
of csoft fields, the csoft function 

Sc&ri = ^5](0|4l4(0)^-n-p)y^(0)|0), 

C 

(13) 

where again all color indices are traced implicity. We 
will see explicitly that both gi/p(£, g) and S c (£, g) indi¬ 
vidually contain rapidity divergences which cancel in the 
total PDF in Eq. (fill) . The appearance of the csoft func¬ 
tion in Eq. (fl3l) is the only deviation with respect to the 
case 1 — x ~ Aqcd/Q, where instead the analogue matrix 
element with soft fields appears, which reads 

S(t, g) = ^<0|4S„(0) S(£-n- p) Sj5 fi (0)|0). (14) 

Note that the structure of the Wilson lines in S c {t, aO and 
S(£, /i) is identical (Xn Sn, only related by a 


In Refs. [T 9 L [20l| these were obtained by resolving the collinear 
sector. Integrating out offshell fluctuations of virtuality ~ Q 2 ( 1— 
x) generate Vn (in analogy to the hard collinear Wilson line Wfi) 
and the Bauer-Pirjol-Stewart (BPS) field redefinition leads to the 
emergence of Xn (in analogy to the ultrasoft Wilson line Yn). 


7 Here the convolution is in fact spurious due to the overall delta- 
distribution in the collinear function g q /p in Eq. (j 12 1) , such that 
the PDF 4> q /p could be also written as a simple product. How¬ 
ever, the form in Eq. m will be more convenient for discussing 
explicit results since the two functions S c and g q /p have the 
same dimension. 
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common boost) and the lack of additional relevant softer 
modes in SCET II_)_ implies that the Lagrangian in the 
csoft sector can be replaced by the full QCD Lagrangian. 
Thus the interactions for csoft and soft modes are equiv¬ 
alent. So S c (£, /j,) and S(£, /i) give the same result (which 
we will demonstrate explicitly in Sec. IV All and we will 
not distinguish them any more in the following discus¬ 
sion of the factorization theorems. Since this concerns 
the only potential difference for the two scaling hierar¬ 
chies 1 — x Aqcd IQ and 1 — x ~ Aqcd /Q, we obtain 

that the factorization theorem in the complete endpoint 
region is always the same for any 1 — x > Aqcd IQ- 
We emphasize that in contrast to the OPE regime 
1 — x ~ 0(1) the endpoint PDF does not encode only 
collinear initial state radiation, but also (c)soft inter¬ 
ference effects between initial and final state radiation. 
Therefore, the PDF at the endpoint may not be inter¬ 
preted only as a description of the momentum distribu¬ 
tion inside the proton before the hard interaction. 

The full factorization theorem reads (to all orders in 
a s and at leading order in 1 — x) 

Fi(x,q) = ^f 2 (x,q)= 

i=q,q 

x J ds J (nj) (^/p ^1 -X- ^2,/^ , ( 15 ) 


where the superscript (rif) indicates the number of active 
quark flavors relevant for the RG evolution of all renor¬ 
malized structures including in particular also the strong 
coupling constant. Here the hard function H ^ nf ) ( Q , q) is 
the square of the matching coefficient between the SCET 
and the QCD currents C^g) (Q, fi) in Eq. ([5|>. while the 

jet function j( n f\s,/j,) and the PDF <j>^Jp(£,n) are de¬ 
fined in Eq. 0 and Eq. m respectively. Note that the 
hadronic tensor becomes transverse in the limit x —> 1, 
such that Fp(x,Q) = 0 and the Callan-Gross relation 
F 2 (x,Q) = 2xFi(x, Q) is satisfied to all orders in a s . 

The massless fixed-order hard and jet functions, 
H( n O(Q^ H ) and (s, hj), are known up to O(a^) 
and O(a^), respectively, and the anomalous dimensions 
are known up to O(a^). Explicit expressions can be 
found e.g. in Ref. |13l | . For the hard function we write 


H^\Q,fi) = 1 + H^’ 1 \Q,n) + 


H, 




+ ffg' 2> (Q,M) + n,frg'" , W,rt +OK), (16) 


r(n/,2) 


where H^’ 2 \ H^’ 2 ' 1 and f’ 2 ’* denote the 

contributions at S(a s ), 0(a?Cp), O(ci^CfCa) and 
0(a^CFTp), respectively. We use an analogous notation 
for all other perturbative expressions throughout this pa¬ 
per. The additional dependence on a finite quark mass 
will be indicated in the arguments. 

The factorization theorem of Eq. is written with 
all its components at the common renormalization scale 


/q which can be chosen independently from the respec¬ 
tive characteristic scales q h ~ Q for the hard function, 
A ij ~ Qy/1 — x for the jet function and ag ~ Aqcd for 
the PDF. Since the choice of fi necessarily differs widely 
from at least two of the characteristic scales, it is manda¬ 
tory to sum large logarithmic terms using the RG equa¬ 
tions. This is achieved by writing each component of the 
factorization theorem as a function that is defined at the 
respective characteristic scale /i//, AG or AG supplemented 
by a RG evolution factor that sums the logarithms be¬ 
tween the characteristic scales and the common scale 

H {nf) (Q,n) 


= 

H^(Q, m ) 

U ( h'\Q,Hh,h), 

(17) 

jG/1 

\s,n) 



= 

J da' J {nf) ( s 

- s',hj) u ( ; f) ( S '^^j) , 

(18) 

<£("/) 

(1 - Z,fl) 



= 

J dz' </> ( g)(G 

- Z, /g) U^ f) (l- . 

(19) 


The individual functions at the respective characteristic 
scales fiHj AG and AG> which are free of any large logarith¬ 
mic terms, serve as the initial conditions of the respective 
RG evolution which follows the RG equations 

A^ U { p f) (Q,^H,^) 

= 7 % f \Q^)U { H f \Q,»H,»), (20) 

A^C/j^W, /j,j) 

= j ds'7j n/) (s-s , ,M)C / j n/) (s',A‘,Mj) ! (21) 

A^ ^ n/) (! - G/q ag) 

= J&z'^ f) (z' - z,n)U^ nf \l-z',^,^). (22) 

The superscript (rif) for all components of the factoriza¬ 
tion theorem (including RG factors) is a reminder that a 
renormalization scheme with rif dynamic running quark 
flavors is used, associated to an n/-flavor scheme. For the 
hard and jet functions and the PDF this scheme is im¬ 
plemented through the common MS subtraction scheme 
for all corrections coming from nj quarks. Subsequently 
this scheme implies that all these quarks enter the RG 
equations via a global nj-dependence. We recall that the 
anomalous dimensions can be determined from the coun¬ 
terterm factors Z that arise in the renormalization pro¬ 
cedure of the individual functions. For massless quarks 


Note the convention concerning the ordering of the arguments 
of the evolution factors for the hard function in comparison to 
the jet and PDF functions. 
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FIG. 3. Exemplary diagrams for secondary massive quark 
production in DIS at O^oICfTf). 


they are defined in the MS scheme. For example, for the 
PDF one has 

dz ' Z ^( z ' - AA 4 ) A^ ~ z '. A 4 ) • 

(23) 

In Eq. (ITSl) the choice of (j, is arbitrary, and the depen¬ 
dence on /j, cancels exactly working to any given order in 
perturbation theory. The fact that any other choice for 
/i can be implemented leads to a consistency relation be¬ 
tween the renormalization group factors, which reads [l3| 

Q 2 U { ?\Q,no,n) Uj nf) (Q 2 (l - 0 ) 

= ( 24 ) 

and a corresponding relation for the anomalous dimen¬ 
sions. 

Accounting for massive quarks the factorization theo¬ 
rem m stays valid with some modifications. This con¬ 
cerns the mass dependence of the hard and jet functions, 
which will be indicated in the arguments, and a modified 
RG evolution with an adapted flavor number according 
to the hierarchy between the renormalization scale n and 
the mass scale /i m ~ m, as described below. The con¬ 
sistency relation (l24l) remains intact since the UV diver¬ 
gences are mass independent. However, additional con¬ 
sistency relations emerge between threshold corrections 
arising when massive quark modes are integrated out. 
In the following sections we will discuss these points in 
detail. 

An important feature of the factorization theorem in 
Eq. (1151) is that there are no flavor mixing terms between 
quarks and gluons in any of the EFT contributions in the 
hard current matching, the jet function, the PDF and 
their evolution factors. One can easily explain this using 
the possible interactions in SCET at leading order in the 
power counting parameter. For example, for the PDF 
evolution one can make the following argument: Flavor 
mixing requires the splitting of an initial state collinear 
quark or gluon into two partons. To stay in the endpoint 
regime one of the final partons has to carry the longi¬ 
tudinal momentum fraction £ > x s=s 1, i.e. almost the 
whole momentum, implying that the remaining parton is 
(ultra)soft. The only available interaction of this kind at 
leading order in SCET is the emission of an (ultra)soft 
gluon from a collinear quark, whereas the splitting of a 
collinear gluon into a collinear and (ultra)soft quark is 


suppressed by 0(1 — x). This means that the parton ex¬ 
tracted out of the PDF at the low scale ~ Aqcd is also 
the one interacting with the hard photon and entering 
the final state jet, and thus cannot be a gluon. Since we 
assume m Aqcd, so that the heavy quarks are not 
produced nonperturbatively out of the proton, this has 
the consequence that massive quarks enter the EFT com¬ 
ponents of the factorization theorem only via secondary 
corrections, i.e. via contributions which are initiated by 
massless quarks and where massive quarks are produced 
through the radiation of virtual gluons that split into a 
massive quark-antiquark pair, see Fig. [3] We mention 
that in the full QCD current there are also flavor mixing 
corrections with massive quarks which in general start 
contributing also at O(al) like the secondary corrections. 
Since these types of corrections do not have a correspond¬ 
ing EFT counterpart, they can be easily included in the 
hard matching coefficient, and we will not consider them 
specifically in our discussion. In fact, due to Furry’s the¬ 
orem, for the case of the electromagnetic vector current 
these effects do not show up at O(o^CfTf) relevant for 
N 3 LL resummation, which is the order where we give 
explicit results. 


III. MASSIVE QUARK CORRECTIONS FOR 
ALL COMPONENTS 


In this section we summarize all ingredients of the 
VFNS for secondary massive quark effects for the most 
singular terms in inclusive DIS in the endpoint region 
x —> 1. We consider a setup with ni massless flavors 
and one massive quark species with mass m 3> Aqcd, 
which we want to incorporate into the factorization the¬ 
orem of Eq. (fl5l) FI This can be easily generalized to 
the case of several massive quark flavors with different 
masses appearing in practical considerations where the 
masses of both the charm and bottom quarks may be rel¬ 
evant. We do not impose any restriction concerning the 
relation of the mass m to any of the hard, jet, or PDF 
scales. The massive quark flavor will never be integrated 
out (in a strict EFT sense) and is thus contained in the 
full QCD description (relevant for the matching compu¬ 
tation for the hard function) as well as in the collinear 
and soft/usoft sectors of SCET. Compared to Eq. (fl5l) 
the factorization theorem for the form factors Fj^ has 


9 We remark that we will not consider the possibility of having an 
intrinsic charm contribution with m ~ Aqcd- In this case the 
mass effects in the perturbative corrections are anyway power 
suppressed. 





an additional dependence on the mass to: 

Fi(x, Q, to ) = ^F 2 (x,Q,m) 

= % H{nf ' ) (Q, m >fJ') j ds J (n/) (s,TO,pi) 

i=q,q 

X ^i/P ~ X ~ ^2 ’ TO ’ ' ( 25 ) 

The factorization theorem has an additional residual de¬ 
pendence on the flavor threshold matching scale fx m . The 
scale fx m is chosen close to the quark mass to, fx m ~ to, 
but otherwise it is arbitrary and represents the scale at 
which we switch between the (m) and the (m + 1) run¬ 
ning flavor schemes. To be specific, we employ for all 
components of the factorization scheme depending on the 
relation of the common scale /i with respect to fi m : 

• the (rif) = (ni) flavor scheme for fx < fx m and 

• the (n,f) = (ni + 1) flavor scheme for fx > fx m . 

These two flavor schemes are implemented independently 
for each of the components in Eq. (1251) by either using an 
OS or the MS subtraction prescription for UV-divergent 
secondary massive quark loop corrections. 

In addition, there are flavor threshold matching condi¬ 
tions arising from the difference of the subtraction pre¬ 
scriptions of the two schemes whenever the RG evo¬ 
lution requires a transition through fx m . The overall 
RG-invariance of the factorization theorem concerning 
changes of the renormalization scales and the consistency 
of properly employing OS and MS subtraction for the 
secondary massive quark corrections ensure that the fac¬ 
torization theorem is continuous at the RG transition 
through the threshold scale \x m . 


A. VFNS for the Hard Function 

The hard function H^ nf ^ (Q, to, h) appearing in the fac¬ 
torization theorem of Eq. (l25l) is in the (nf) = (ni + 1) 
flavor scheme if fx is above /x m and in the (n/) = (rq) 
flavor scheme if fx is below /r m : 

f H( ni+1 \Q, to, fx) for fx > ix m , 
HM(Q,m,fx)={ 

[HM(Q,m >f x) for /x < fx m . 

(26) 

The common renormalization scale /r is in general dif¬ 
ferent from the characteristic scale /iff ~ Q of the hard 
function, so we specify the hard function at the common 
scale fx by the hard function at the scale /iff, which is 
free of any large logarithmic terms, multiplying a RG- 
evolution factor that resums the logarithms between the 
scales /iff and /i. The hard function at the scale /iff 
serves as the initial condition of this RG evolution, and 


the flavor scheme that is employed for the initial condi¬ 
tion depends again on the relation of /iff to the flavor 
matching scale ix m : 

[ H( ni+1 \Q,m,Lx H ) for/iff >/t m , 
H^ f \Q,m,ix H )= l 

{ H( n, \Q,m,n H ) for ^ ^ Mm 5 

(27) 


where 

H^Xq, to, /iff) = H^XQ, /iff) + 21 ^ 02 ) (Q) m); 

(28) 

H {n ‘ +1) (Q, to, /iff) = tfOu+i) (Q, /iff) + 2F^ +1 ’ 2) (Q, to ). 

(29) 


The functions H^ ni XQ, I^h) and H^ ni+1 XQ, /iff) are the 
hard functions for massless quarks in the ( ni ) and (?ij-|-l) 
flavor scheme respectively in the notation of Eq. m- 
The term Fm 1 ’ 2 ' 1 represents the massive quark loop con¬ 
tribution to the QCD current form factor with OS sub¬ 
traction in the (m) scheme, see the first diagram in Fig. [3] 
In fact, for the case of DIS in the endpoint region all 
mass dependent corrections at O^oICfTf) can be di¬ 
rectly inferred from the matching calculations carried out 
for thrust in Ref. @ due to the fact that the hard coeffi¬ 
cients are the same up to an analytic continuation from 
the timelike to the spacelike process. We write the result 
as (ai' ll) = ai ni) (/rff)) 


F£" 2 XQ,rn) = 


>0\ 2 


) CfT f 


167r 2 


(30) 


(2\ 

with the function fm (Qi m ) give 11 by 



where 

r = Vl - 4 to 2 , (32) 


and to = m/Q. Due to the imposed OS renormaliza¬ 
tion condition this correction decouples for large masses, 
i.e. Fm l ' 2 XQi m ) —t 0 for to = m/Q oo. For small 
masses the function ( Q , to) reads 


/i 2) (Q, ™) c j In 3 (to 2 ) + y In 2 (to 2 ) 


+ 


/ 530 4tt 2 

V~27 + ~ 9 ~ 


2 \ .. 2 , 3355 

1ii(to 2 ) H- 

V 2 81 


38t r 2 

~2T~ 



(33) 


which exhibits mass singularities and is therefore 
not suitable in the small mass regime. The term 
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^ +1 ’ 2) (Q,H represents the corrections due to the 
nonvanishing mass of the heavy quark in the (rq + 1)- 
flavor result and can be cast into the simple form 


p(n, + l,2) 
^ Am 


(Q,m) 


(ai ni+1) ) 2 C F T F 

I6n 2 



m—> 0-1 


(34) 


which can be read off explicitly from Eqs. (1311) and (1531) . 
Due to the fact that F^ +1,2 ^(Q,m) —»• 0 for m = 
m/Q ^ 0 the massless limit is recovered in the hard func¬ 
tion for a vanishing quark mass. The form of Eq. (1341) is 
a direct consequence of the fact that the leading IR de¬ 
pendence has to cancel in the SCET matching (in MS). 

The RG evolution from /j, F to the common scale /r in 
the VFNS proceeds in the (rq + 1) flavor scheme as long 
as the scale is above and in the (rq) flavor scheme if 
the scale is below /i m , according to Eq. COl) . Finally, if 
the RG evolution crosses the flavor matching scale fi mi 
one has to account for a threshold correction factor which 
we call if RG evolution crosses from the (rq) flavor 
scheme to the (rq + 1) flavor scheme and Ai ^ if RG evo¬ 
lution crosses from the (rq + 1) flavor scheme to the (rq) 
flavor scheme. They are the inverse of each other since 
they are just the ratios of the hard function of Eq. (EH 
in the two flavor schemes: 


to, n m ) 


M H (Q,m, lim) 


H( ni+1 \Q,m, n m ) 
H( n ‘'>(Q,m,)j, m ) 
ff (ni) (Q,m,/r m ) 
77(«<+ 1 )(Q,m,/z m ) ' 


(35) 

(36) 


Here the ratios should be expanded with a common 
choice of either ai ni \ii m ) or (/r m ). The thresh¬ 

old correction factor for the hard function at fixed 0(a 2 ) 
reads ( L m = In (m 2 / fi 2 )) 




FO 


a 2 C F T F 
167r 2 


x 


r 2 80 r 
Lm + —L m + 


3 
260 


+ 


47+ 


Lm + 


224 

~27 

875 


In 


m 

Q 2 

107r 2 



where a s can be either written in the rq or rq +1 scheme 
and = —M~ I } 2 \ The corrections at fixed 0(a s ) 

are zero, i.e. fj, m ) = 1 + A4^ (2) +0(ofy. Note 

that the threshold corrections A4^ involve the logarithm 
ln(m 2 /Q 2 ). It is directly related to the rapidity diver¬ 
gences arising in the computation of the SCET current 
with collinear and soft fluctuations tied to the mass shell 
of the secondary massive quarks, see Ref. jhj. In the loga¬ 
rithmic counting 0(a s ln(m 2 /Q 2 )) ~ 0(1) one therefore 
has to include also the terms at 0(o? s \n(m 2 /Q 2 )) and 
0(a^ ln 2 (m 2 /Q 2 )) for a computation of Mg at preci¬ 
sion 0(a 2 ). In Sec. [V] we will show how to obtain these 
terms via consistency from the PDF threshold correction, 


for which we will also demonstrate explicitly the expo¬ 
nentiation property of the rapidity logarithms. 


B. VFNS for the Jet Function 


The VFNS for the jet function can be set up in a 
way analogous to the hard function. The jet function 
j( n O(s , m, /r) is in the ( rif) = (ni + 1) flavor scheme if /i 
is above and in the (ny) = (rq) flavor scheme if /j, is 
below ii m : 

f j( ni+1 )(s, m, n) for n>nm, 
j( nf \ s ,m,n) = < 

[ (s, m, n) for n < n m ■ 

(38) 

The common renormalization scale /r is in general differ¬ 
ent from the characteristic scale [ij ~ Q\J 1 — x of the jet 
function, so we specify the jet function at the common 
scale /r by the jet function at the scale /rj, which is free 
of any large logarithmic terms, convoluted with an RG- 
evolution factor that resums the logarithms between the 
scales /ij and p,. The jet function at the scale p,j serves as 
the initial condition of this RG evolution, and the flavor 
scheme that is employed for the initial condition depends 
again on the relation of fj,j to the flavor matching scale 
f^m • 

[ J (n!+1) (s,TO,/7j) for [lj > Hm , 

^ n/ ( s , m, hj) = < 

[ J (7li) (s,m,/xj) for [ij < Hm , 

(39) 


where 

J M (s,m,fi) =J ( ' ni \s,n) ( 4 °) 

J ( ” i+1) (s, m, n) = J ( ” !+1) ( s , n) + (s, m,/x) 

+ J ( Zt* 2) M- ( 4 !) 


All mass dependent corrections at 0(a 2 C F T F ) can be di¬ 
rectly inferred from the results computed in Ref. [@J due 
to the fact that the thrust jet function is decomposed 
out of two hemisphere jet functions each of which are 
the same as the one in DIS. The terms ./^"('^(.s, m) and 
Jm ‘ reai’ 2) ( s , m ) in Eq. (EM contribute only when the jet 
invariant mass is above the threshold 4 to 2 and thus cor¬ 
respond to real production of the massive quarks. They 
are given by 


T (n/,2) 1 , (a { r f) ) 2 C F T F l of a 

J mieal(+ TO ) = - T7r~ 2 - - 9( S - Am ) 


32 

— - Ll2 

3 


6-1 
1 + 6 


167T 2 

16 , 
+ y ln 


- - In" 


1-6 

lTfe 


1 — 6 " 
4 

- 6 4 — 6 2 + 

2 18 


In 

In 


1-6 
T+b 
1 - 6 
1 + 6 


5 ,, 241 , 8 tt 2 

- b 6 -\ -6- 

27 9 9 


( 42 ) 
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for both n/ = rq and rif = ni + 1 with 



Note that is zero at the threshold s = 4m 2 

and that it decouples for m —> oo automatically due to 
the threshold 0-function, as required by the OS prescrip¬ 
tion. The expression for J ^ + c ^ (s,m, n) contains only 
distributions and corresponds to collinear massive vir¬ 
tual corrections (including soft-bin subtractions) as well 
as terms related to the subtraction of the massless quark 
result already contained in J^ ni+1 \s, n). Its renormal¬ 
ized expression reads (s = s//z 2 ) 


V 2j A£d£t( s ’ m ’V) = 

58 r 2 / 718 8 tt 2 

+ ^L 2 + I - 


(ai ni+ 1 ) ) 2 C F T F 


\ 27 


9 


L r , 


167T 2 
4325 


58t r 2 


_ r 3 

9 m 


e( s ) 


J J+ L 
9(s) Ins 


_8 2 _ 116 r 
3 m 9 m 27 


81 

718 8 tt 2 


_32, 

27 3 


- 1 + L 


16 116 
—L m + — 


9(s )In 


2 —' 


(44) 


The jet function in the (rq + 1) scheme reaches the mass¬ 
less limit, i.e. j(”' +1 )(s,m, /z) —*■ j( n ‘ +1 )( S) ^) for m —> 0. 

The RG evolution from /xj to the common scale /i in 
the VFNS proceeds in the (rq + 1) flavor scheme if the 
scale is above p, m and in the (rq) flavor scheme if the scale 
is below [i m , according to Eq. (Ell) . Finally, if the RG 
evolution crosses the flavor matching scale /z m , one has 
to account for a threshold correction factor which we call 
Aij if RG evolution crosses from the (rq) flavor scheme 
to the (rq +1) flavor scheme and Mj if RG evolution 
crosses from the (rq + 1) flavor scheme to the (rq) flavor 
scheme. They are the inverse of each other since they 
are just convolutions of the jet function of Eq. (1391) in 
the two flavor schemes: 


Afj'(s,TO,/r m ) (45) 

= (s',m,/z m ), 

Mj (s, m, Hm) (46) 

= (s',m,/im). 


The threshold correction factor for the jet function at 
fixed 0(a 2 ) reads 




58,2 
9 Lm 


9(s) 


466 


-I + L 


Ll + 


FO 

47T 2 

80 


a 2 s C F T F 
16-7T 2 
1531 


Lm. 


S(s) 

107T 2 


__ t3 

9 Lm 


54 


27 ' 9 


224 

9 m + ¥ 


(47) 


where a s can be either written in the (rq) 
or (rq + 1) scheme and A4j^ = 

The corrections at fixed 0(a s ) are zero, i.e. 
M j (s,m, /im) = S(s) -I- A4 + O(al). Note that 
implicitly contain a logarithm ~ ln(m 2 /s) that 
becomes large for m -C y/s ~ Q\J 1 — x or m y/s. Its 
presence becomes more manifest when using the natural 
scaling variable s = s/v 2 ~ 0(1) with vj ~ Qy/1 — x 
instead of s = s/ 


v jM J ( - 2 \s,m,iJ irn ,uj) = 


58 2 / 466 Rr 2 

" r m ” V"27 + ~~9~~ 


+ 


0(S) 


In 


-i + 



In the logarithmic counting 0(a s ln(rrr 2 /s)) ~ 0(1) one 
therefore has to include also the terms at 0(a^ In (m 2 /s)) 
and 0(aj In 2 (m 2 /s)) for a computation of M J at 0(a 2 ). 
The corresponding results will be discussed in Sec. |V] 
where we show how to obtain these terms via consistency 
from the PDF threshold correction. 


C. VFNS for the PDF 

Analogous to the case of the hard and the jet func¬ 
tion one can set up the VFNS for the PDF. The PDF 
cj)( n f)(l — z,fj) is in the (n/) = (rq + 1) flavor scheme if 
/z is above /z m and in the (n/) = (rq) flavor scheme if /z 
is below n m : 


\ 0("‘ +1 )(l -z,m,n) for /z>/z m , 
4> {nf) ( 1 - z,m,n) = < 

[ 0( nj ) (1 - z,fj,) for fj, < . 

(49) 

Note that in the (rq) scheme the dependence on the quark 
mass vanishes for m Aqcd- The common renormaliza¬ 
tion scale /j, is in general different from the characteristic 
scale /fo ~ Aqcd of the PDF, so we specify the PDF at 
the common scale /x by the PDF at the scale /z^, con¬ 
voluted with an RG-evolution factor. The PDF at the 
scale serves as an initial condition for the RG evolu¬ 
tion and since we consider effects of heavy quarks with 
a mass m Aqcd, the mass scale is always above the 
scale of the PDF, i.e. /z m > /z^, so that the PDF at /z^ is 
always in the (rq) flavor scheme. This is independent of 
the scaling of (1 — x) with respect to Aqcd/Q- 

The RG evolution from /x^ to the common scale /z in 
the VFNS proceeds in the (rq + 1) flavor scheme if the 
scale is above /z m and in the (rq) flavor scheme if the 
scale is below /z m , according to Eq. (1221) . Finally, if the 
RG evolution crosses the flavor matching scale /z m , one 
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has to account for a threshold correction factor AlJ, 


z,m,/j,m) (50) 

= Jdz' (f>( ni+1 \z' — z, m, fj, m ) (1 -z',m,n m ). 


i = Q( 1 - z) and ~ Q(1 - z ), 


^0 .. + (2) { £ ^ 
— M 4, 


Q 

X J<5(/) 


8 2 80 
3 Lm + Y Ln 


FO 

224 

"27 


o^CfTf 
(47r) 2 


In 


Q 


(2 8tt 2 5 

r 73 

20t r 2 

8 1 

U + - y 

Lm + 18 + 

27 

~ 3 C 3 


'o(i) 


i 

+ L 




80 




224 

77 


(52) 


Since we always assume (jL m > /z^, only the transition 
from the (m) to the (ni + 1 ) flavor scheme is relevant for 
the PDF. The results for A4j(l — z,m,fi) in the end¬ 
point region can be easily obtained from the well-known 
PDF threshold factor in the OPE region calculated in 
Ref. j{| and for convenience also given in Eq. (IA11I) of 
Appendix [A] by expanding for z —t 1, which yields 


The large logarithms ln(l — z) arise from rapidity di¬ 
vergences in the collinear PDF function g i / P and the 
soft function S (or csoft function S c ) and can not be 
resummed in the usual RG evolution. In Sec. m it will 
be shown how these logarithms can be resummed using 
rapidity RG methods as in Refs. mm enabling us to 
evaluate the matching coefficient AlJ at N 3 LL order. 

IV. PRACTICAL IMPLEMENTATION AND 
CONSISTENCY RELATIONS 


- A/{ 1 (2) ( 1 - Z,m,g,m) 


FO 


u 2 s C f T f 

(47t) 2 


6(1-z) 2 Ll 


+ G 


8 

3 ' 9 


73 20t r 2 8 

m + 18 + 27 ~ ~i 


+ 


6(1-z) 


1 - 2 


8 2 80 r 

—L -f- — L n 
3 m 9 


224 

77 


( 51 ) 


In this section we specify explicitly the RG properties 
of the individual matrix elements in the factorization the¬ 
orem of Eq. (1251) in the presence of massive quark correc¬ 
tions for all hierarchies between the mass scale and the 
kinematic scales. All of the factorization theorems are 
valid up to power corrections of 0(1 —x), independent of 
the hierarchy between the mass and the kinematic scales, 
since the change between the OS and MS renormalization 
prescriptions does not generate any power corrections in¬ 
volving the mass. Furthermore, we investigate the con¬ 
ceptual implications of using different final renormaliza¬ 
tion scales which all of the matrix elements are jointly 
evolved to, which leads to consistency relations among 
the threshold corrections. 


where a s can be either written in the (m) or (ni + 1 ) 
scheme. The corrections at fixed 0(a s ) are zero, 
i.e. M^(l-z,m,ii m ) = 6(1-z) + +0(a 3 s ). In 

Sec. |V] we will also compute this result directly from the 
definition in Eq. m for the PDF in the endpoint re¬ 
gion. Ait contains a large logarithm ~ ln(l — z) that is 
manifest when rescaling the plus-distribution in terms of 
the normalized soft momentum variable t = tj with 


A. Explicit factorization theorems with massive 
quarks 

We apply the prescriptions given in Sec. IIIII and use 
first for definiteness a common renormalization scale /i < 
g m . This implies that the matrix elements and couplings 
are renormalized in the ( ni ) scheme at the common /i. 
If the RG evolution from the natural scale of the matrix 
element /j,; (for i = H, J) to the final scale g, crosses the 
scale [u,rn , the scheme is changed leading to the threshold 
correction M~ and the number of active flavors in the 
evolution changes from ni + 1 to nj. Thus the hard and 
jet functions can be written as 
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H^ ni+l \Q,m,ii H )U^ l+l \Q,ii H ^ m )M H (Q,m,ii m )U ( ^ l) (Q,ii m ,iJi) for /j, H > /i m , 
H^(Q,m,fi) = { (53) 

H( ni \Q,m,ii H )U£ l \Q,iJ, H ,ii) for /x ff < /x m , 


/ d s'f d s” f d s'" J( ni \s - s', m,nj)U^ ni \s' - s", n m , fij) 
J {ni \ S ,m,ll) = l X Mj (s" - S'" ll m ) (s'" - S", [I, /Ira) 

J ds' J (ni) (s - s', to, /xj) Uj ni \s', n, iij) 


for hj > Hm, ( 54 ) 
for flj < Hm ■ 


Note that the PDF is never evolved in the (?x; + 1) flavor scheme for /x < /x m because we always assume to ~ /x m > 
/i 0 ~ Aqcd, so that Eq. (1191) holds with n/ = rq. The explicit description of the complete factorization theorem with 
all evolution factors written out thus adopts three different forms depending on the hierarchy between /x m on the one 
hand and hh and fij on the other. For simplicity we set here fi = n<j>, so that the RG factor for the PDF can be 
dropped. For ii m > hh we get 

F\(x,Q,m) = ^2 yif (ni) (<3 ,iti,iih)u2 i) (Q,IJ-h,ii^) jdsjds' J {ni) (s',m,nj) E/j ni) (s - s'^/xj) 


i=q,q 




1 — X — 


Q 2 


(55) 


This factorization theorem covers in particular the region, where the massive quark decouples and therefore for all 
renormalizable quantities the OS scheme is used for the secondary massive quark effects. In this regime mass effects 
in the jet function are power suppressed by O ^ ~ ^(1 — x ) and might be dropped due to their small sizeF^I 

For hh > l^m > one gets 


F?(x,Q,m) = J2 ^H^ +1 \Q,rn,ii H )U^ l+1) {Q,ii H ^ m )M- H {Q,m,ii rn )U^\Q^ m , H ) 

i=q,q 

x J ds J ds ' • ( 56 ) 

Here we use both the MS and OS renormalization prescriptions for the secondary massive quark effects in the evolution 
of the hard function. In particular the MS scheme allows us to reach the massless limit for the fixed-order hard 
coefficient H( ni+1 \Q,m, iih)- Since /x m > fij the jet function is still always evolved in the (m) scheme. It is easy to 
see that for /x m = /x h the two factorization theorems in Eqs. (l55l) and (15B1) agree due to the matching relation (1561) . 
So there is a (perturbatively) continuous transition between the two scaling hierarchies described by and Fj IT . 
Finally, for /xj > fi m > /x^ the explicit factorization theorem reads 

2 

*?“(*, Q,m)=Y,T H(n,+1) (Q, m, rn) U { h 1+1 \q, m, to, /x m ) U { ” l) (Q, /* m , h) 

i=q,q 

x J ds J ds' J ds' 1 J ds'” J (ni+1) (s'", to, /xj) ( s " - s'”, /x m , /xj) Mj (s' - s”, m, /x m ) 


X (s - s', ll,,/,, Hm) (jtyp ^ 


1 — x — 


Q 2 


(57) 


10 Since mass effects in the jet sector appear only in the real radi¬ 
ation correction ./((((.^ containing a kinematic threshold, they 
can in practice anyway not contribute. 


Now in addition both renormalization prescriptions 
are also used for the jet function allowing us to 
reach the massless limit for the fixed-order structure 
j( n i+i) (g, to, yixj). In this regime mass corrections in 

the hard function are power suppressed by O (fp) < 
0(1 — x ) and taking the massless limit might be suit- 
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able. Again, it is easy to see that for /i m = fjj the two 
factorization theorems in Eqs. and (1571) agree due to 
the matching relation (Idfil) . So there is a (perturbatively) 
continuous transition between the two scaling hierarchies 
described by Ff 1 and F™. 

B. Consistency conditions 

The equivalence of the factorization theorem for differ¬ 
ent choices of the common renormalization scale /r con¬ 
cerning physical predictions leads to statements about 
the intrinsic relations between its components. On the 
one hand, they imply the well-known consistency condi¬ 
tions between the RG evolution factors Uh, Uj and U 
see Eq. (1771) . On the other hand, in the context of the RG 
evolution crossing a massive quark threshold they also 
imply a consistency relation between the threshold fac¬ 


tors for the hard, jet and parton distribution functions, 
M.%, Af j and VU J. Apart from providing consistency 
checks of theoretical calculations, these relations have 
also computational power, as they can be used to calcu¬ 
late properties of independent gauge-invariant field theo¬ 
retic objects once it has become clear that they represent 
building blocks of a factorization theorem. Hereby, one 
of the most interesting aspects is that the various build¬ 
ing blocks can appear in different factorization theorems, 
and one may gain insights into the mass-singularities of 
apparently unrelated quantities. 

In the previous subsection we have discussed the renor¬ 
malization of the matrix elements for /r < /j, m . An equiv¬ 
alent choice would have been fj, > /x m , where all renor¬ 
malized quantities are evaluated in the (nj + 1) scheme at 
the common scale /i. Here the RG evolution of the hard 
and jet functions and the PDF reads 


H^ +1 \Q im , (jl) 


# (rii+ 1 ) (Q, m> ^ H ) u%' +1 ) {Q,h h , fi) 

-ff (,l!) (Q, m, ij, h ) U£'\Q, uh, lira) Mff(Q, to, /j m ) uj?‘ +1) (Q, n m , /r) 


for /iff > Hm, 

(58) 

for /iff < Hm , 


( f ds'J (ni+1) (s - s',m,/j,j) U^ ni+1) (s',fj,,/ij) 

J(n ‘ +1) (s ’ = ) f d s' / ds" / ds" f d s"'JM{8 - s', m, fij) Uj ni \s' - s", /x m , /ij) 

{ x Mj(s" - s'", m, n m ) Uj ll+1 \s" 1 - s",/i,/i m ) 

0 (n ' +1 )( 1 - z, m, fj.) = J dz'dz"dz'"</> (n,) (z' - z,/*«/,) U^ ni \z" - z', /x m , /i</>) 

X M J (z'" - z", m, Hm) U^ l + 1) (1 - z'", H, Hm) ■ 


for UJ > Hm , 


for UJ < Hm, 


(59) 


(60) 


Here the threshold factor (1 — z, to, /i m ) arises since the RG evolution of the PDF </) necessarily crosses the massive 
quark threshold. 

In Fig. [I] we show an illustration of the two equivalent choices for fij > /i m , i.e. we display the situations where 
the common renormalization scale /i lies (a) between the mass and the PDF scales or (b) between the jet and the 
mass scales. To discuss the implications let us consider the complete factorization theorem for this specific hierarchy 
in case (b) with /i = fjj set for simplicity, 


Fl u (x,Q,m) = ^ ^-H ( ” i+1) (Q,m,//ff)I/^” i+1) (Q,/iff,/xj) J ds (s,m, uj) J dz U^ l+1) (z - x, uj, fi m ) 


i=q,q 


x J Az ' j d z"M^{z' - z,m,Hm ) U^ ni \z" - z', fj, m , (l - z" 


Q 2 ’ 




(61) 


The equivalence of the factorization theorems in Eqs. m 
and (l6ll) implies, besides the relation between the evolu¬ 
tion factors and anomalous dimensions shown in Eq. (1241) 
for rif = ni and ni +1 , also a relation between the thresh¬ 
old correction factors, 

M^( 1 - z,m,n) =Q 2 Mx(Q,m : n)Mj(Q 2 ( 1 - z),m,/x) 

(62) 


or equivalently 

d(l - z) = Q 2 Mff(Q, m, n) (63) 

x J dz' Mj(Q 2 (l - z'),m,/i)Mj(z' - z,m,n). 

These relations imply in particular that the rapidity loga¬ 
rithms (and singularities) that arise in the hard, collinear 
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(a) (b) 

FIG. 4. Illustration of the different RG setups for the hierarchy /xj > /x m > /x^, leading to the consistency relations mentioned 
in the text. We display the cases where the common renormalization scale /x satisfies (a) /x m > /x > /x,/, and (b) gj > /x > /x m . 


and soft sectors are intrinsically related to each other. 
We can explicitly check that the consistency relation is 
satisfied at O(o%) in the fixed-order expansion. Inserting 
Eqs. ([37jl . (1471) and ACT) confirms Eq. (1521) . We emphasize 
that for Eq. (TSTT) to be satisfied for arbitrary masses the 
coefficients of the rapidity logarithms ln(m 2 /Q 2 ) in M^, 
ln(Q 2 (l — x)/m 2 ) in M ± and ln(l — x) in need to 
be equivalent. 


V. CALCULATION OF THE PDF THRESHOLD 
CORRECTION 

As already indicated in Eq. (15U1) . the PDF threshold 
correction A4^(l — z,m,/x) is given by the ratio of the 
PDFs in the (rq + 1) and (rq) schemes, i.e. 

1 - z, m, n m ) 

= j A z'^\l-z',m,, m) {^) (z' - z,m,iu, m ) 

= Jdz'Z { ™ l \l - z',m,n m ) (z'-z,/x m ), 

(64) 

where the second equality arises from the universality of 
the unrenormalized bare PDF. Note that in Eq. (1641) the 
calculation of the PDFs can be performed with partonic 
initial states (i.e. quarks) since the different renormaliza¬ 
tion conditions are not affected by the infrared behavior. 
This can be also seen from the second equality which only 
involves the renormalization factors. 


Since the PDF in the endpoint region is decomposed 
out of a soft function S (or a csoft function S c , depending 
on the scaling of (1 — x) with respect to Aqcd /Q) and a 
collinear function g , see Eq. the analogous relations 
to Eq. (1511) hold also for the corresponding matchin g co - 
efficients Ms and M g , which are related to vhrl 


A4j(l - 2 :,rn,/x m ) = Q J (65) 

X A 4g (Q(l Z ) f, (Ji 1X1 , /Xm, xx) . 

An important technical point is that we encounter rapid¬ 
ity divergences in the calculation of the collinear and soft 
PDF functions which are not associated to the UV or IR 
behavior and are not regularized by dimensional regular¬ 
ization. To regulate these divergences we need an addi¬ 
tional regulator that breaks boost invariance. Here we 
display the corresponding results for individual diagrams 
employing the “^-regulator” [lo|.[TT| for the collinear, soft 
and csoft Wilson lines, i.e. 


W n = ^ exp 

perms 


9 \n ■ V\~ v 
n • V 


^ * A n 7 


( 66 ) 


Sn = ^2 eXP 

perms 


9 |2P 3 |~ t?/2 

n ■ V 


n- A s , 


(67) 


11 For notational simplicity we use here the soft matching coefficient 
A4 s, which is identical to the csoft matching coefficient Ais c • 
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X n = '^2 exp 

perms 


g \n ■ r P\~ r i/ 2 
n ■ V i /-*?/ 2 


n • A 


CS 


( 68 ) 


and similarly for S n and Vn, where due to the boost of 
the csoft modes with respect to the soft modes \ 2 V 3 ,\ is 
replaced by \n ■ V\. In this context the scale v is an aux¬ 
iliary scale to maintain the dimensions of the regulated 
integrals which adopts a similar role as the g scale in di¬ 
mensional regularization. In particular, also the strong 
coupling adopts a ^-scaling proportional to g. 

We follow the method of Refs. |To|, [ll| for setting up 
the rapidity RG evolution. The summation of the rapid¬ 
ity logarithms can be carried out independently after the 
/i-evolution has been settled which is the approach we 
are adopting here. We will show that the decomposition 
in Eq. (1651) provides a way to resum rapidity logarithms 
~ ln(l — x) in terms of a RG evolution in v. A similar fac¬ 
torization in rapidity is used in the hard current match¬ 
ing computation for massive primary quarks of Ref. [23j] 
which also discusses the rapidity RG evolution due to 
secondary massive quark effects in detail. 

For sufficiently inclusive observables dispersion rela¬ 
tions can be used to obtain the results for secondary mas¬ 
sive quark radiation (with mass m) at 0(a 2 C F T F ) from 
the results for “massive gluon” radiation (with mass M) 
at 0(a s ), which allows us to deal with the technically 
simpler one-loop computations for the latter instead of 
performing the two-loop integration directly. The dis¬ 
persion method has been discussed in detail in Ref. @ 
and we refer to Sec. IV A therein for the notations and 
the explicit relations involved. For the following compu¬ 
tations we use always Feynman gauge. 


A. One-loop results for the PDF soft and collinear 
functions with a massive gluon 

For the computation of the PDF threshold correction 
we have to consider both the collinear PDF function 
g q /p(£ : g) defined in Eq. IT2l) and the csoft or soft func¬ 
tions S c (£. ji) and S(£,p) defined in Eqs. (fl3l) and (fT4l) . 
respectively. As already argued in Sec. Ill Bl the csoft and 
soft function are related by a common boost of the Wilson 
lines and are therefore in fact equivalent, which we will 
explicitly show here at the one-loop level. Since we are 
interested in the matching correction related to different 
employed renormalization schemes, we can perform the 
computation with partonic initial states. In this subsec¬ 
tion we consider only the massive gluon contributions. A 
similar calculation has been performed in Ref. ea within 
the context of using the gluon mass as an IR regulator. 

Let us start with the computation of the partonic 
collinear contribution g^ q - Since this is a local matrix 
element, no real radiation diagrams can contribute^ 


= 



FIG. 5. Feynman diagram for the collinear PDF function 
with a massless quark field in the initial state and a massive 
gluon at 0(a 3 ). The symmetric diagram and wave function 
corrections have to be added. 


Therefore the only contribution (besides the wave- 
function renormalization) for massive gluon radiation 
at 0(a s ) is given by the virtual gluon contribution 
of the diagram in Fig. [5] and its symmetric configu¬ 
ration, which we denote together by g n .. For conve¬ 
nience, we use a frame where the perpendicular com¬ 
ponent of the initial onshell quark momentum vanishes, 
i.e. = ( p + ,p~,p ± ) = (Q,0,0). We then obtain 

gn = Aig 2 C F '^vH(£) J ( k+ y+v (69) 

1 1 

[k+k~ -k 2 ± - Qk~ + ie] [k+k~ -k 2 ± -M 2 + ie] ‘ 


After performing the k - and fcj_-integrations we get 


g n= -^ 5 {i) r( 2 -| 

7r \ z 


p, 2 e 7E 

M 2 


2 -§ 


r Q d k+ 

lo 


k + 

1_ ~Q 


(70) 


This contribution is not regularized by dimensional reg¬ 
ularization alone (i.e. for g = 0 ) and there is no collinear 
real radiation contribution to cancel the correspond¬ 
ing rapidity divergence (in contrast to the OPE region 
1 — x ~ 0(1)). Expanding for g —> 0 gives 


a s C F 

9n — 

7r 


s(£) r 2 


d\ //x 2 e 7E \ 2 2 
2 ) M 2 ”") 



(71) 


where H a = 1 + a) + 7 F is the Harmonic number. 

The corresponding soft-bin subtractions gn,0M have to 
be taken into account, which in general yield some ad¬ 
ditional corrections. However, for the 77 -regulator they 
vanish. Including the contribution from the wave func¬ 
tion renormalization, given by 


v(l) a ' C F-nf n d\f^^ 2 -^2(d-2) 
_ ~~Ap~ — 2y v - 


(72) 


need to be soft, a contribution that is excluded from the collinear 
12 To stay in the endpoint region 1-2<1 the emitted gluon would matrix elements by zero-bin subtractions. 
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FIG. 6 . Non-vanishing Feynman diagrams for the computa¬ 
tion of the one-loop massive gluon contributions to the soft 
function S(£, g). The corresponding symmetric configurations 
are implied. 


we obtain in total for the bare partonic collinear function 


= (dn - gn,OM ) - 5(1) 


u s Cf 

47T 


<S(£)r 2 


d\ (n 2 e 7E \ 2 = 

2) \ ) 


+ 4 In 


4JTd_ 1 - 


2(d — 2 ) 


(73) 


This is in agreement with the results given in Refs. (llL 

13- 

Let us turn to the computation of the (partonic) soft 
matrix element S^\ where the corresponding diagrams 
are shown in Fig. [51 The virtual diagrams £ v irt with the 
77 -regulator employed for both Wilson lines read (includ¬ 
ing symmetric configurations) 


S Virt =~Aig 2 C F g 2e v r >5(£) 


X 


d d k | fc+-fc-p _ 

(27r) d [k~ — ie] \k + — ie] [k 2 


(74) 

1 

M 2 + ie] 


A simple way to evaluate this is to first perform the 
/co-integration by contours and afterwards the remain¬ 
ing k$- an d fc 1 - integrations, which yields in agreement 
with Refs. fnlTpi 


-hvirt — 


mr 2 -- 


a s C F 


x { - + 2 In 
.V 


g 2 e' 1E 
M 2 


+ < 75 > 


The real radiation diagrams .SVeai yield 
d d k \k + - 


S rea i = 4 ig 2 C F ~^^ 


(2ir) d [k — ie] [k + — ie] 
x (-2 t ri)6(k 2 - M 2 )5(t~k+). (76) 


After performing the trivial k + - and fcy-integrations this 


reads 

*Sreal = 


a s C F [g 2 e lE Y 


< M2 k 


r(| - 1 ) 1 

Ak - (n k --M 2 y- 2 \t-k 


-\-v 


(77) 


Finally, expanding in 77 after the k -integration yields 
(with 1 = l/v) 




Ot s C F ^ 

7T 


d\ //x 2 e 7E \ 2 = 

2) yJd 2- ) 


(78) 


X 



- + 2 In 
1 





We note that our computation of SVeai differs from 
Ref. 13 which uses the same regularization methods in 
several ways: First, our prescription for the Wilson lines 
in the soft function differs from theirs resulting in a rel¬ 
ative sign in Eq. ( 1751 ) . Second, the result of the phase 
space integrations in Eq. (35) of Ref. [3 does not agree 
with Eq. (1781) . Third, we emphasize that in the compu¬ 
tation of the soft diagrams we do not encounter any non¬ 
vanishing collinear-bin subtractions, in contrast to such 
a statement given there. However, overall these three de¬ 
viations cancel each other giving the same result for the 
total soft real radiation correction in our Eq. ( 1751 ) . 

Summing up all contributions the bare soft function 
reads in terms of £ = £/with v<f, ~ Q( 1 — x) 


V* ^ (bare d } (£, M , (i, v) = V+ (S virt + S rea l) 

a s C F ^(^ d\(g 2 e^\ 2 --* 

4tt ^ 2 )\M 2 ) 

(79) 

in agreement with the result in Ref. 13- 
Next we will also calculate the csoft function S c to show 
that this leads to the same result as for the soft function S 
above. Here only the rapidity regularization prescription 
changes, see Eq. ( 1551 ) compared to Eq. ( 1571 ) . This gives a 
scaleless contribution for the virtual diagram, such that 
S'c.virt = 0 . The real radiation diagrams for the csoft 
function then give 


x 5(1) 


■- -41n( — 

77 W/J 


m 


^c,real — 


a s C F (/i 2 e 7E ) 


2—7 


m 

n T (| - 1) (. 1+r > 


dfc 


' m 2 k 

i 


a s Cp 


( ik- - M 2 ) 
d 


r 2- 


g 2 e JE 

M 2 


2-i 


0 (£) zA 
p+v 


(80) 


which slightly differs from Eq. ( 1771 ) concerning the depen¬ 
dence on the 77 -regulator. Expanded in 77 gives exactly the 
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same result as shown in Eq. d- Finally, we remark that 
at one loop the 77 -regulator for the csoft function acts in 
the same way as the a-regulator suggested in Ref. [24| 
applied to the large light-cone component, 


dfc+ a dfc+ 
k+ (k+) 1 + a ' 


(81) 


Since the a-regulator is boost independent, it gives au¬ 
tomatically the same result for the virtual and real ra¬ 
diation diagrams of the soft and csoft function. For the 
collinear function g n both regulators are anyway identi¬ 
cal, which thus also implies that the 77 -regulator needs to 
yield the same result for the soft and csoft functions. We 
will therefore not distinguish between the soft and the 
csoft functions anymore in the following. 

Expanding Eqs. d and d> for d —> 4 gives the un¬ 
renormalized correction^ ( Lm = In(M 2 //x 2 ), £ = l/v$) 


g^ ie ’ 1 \£,M,Q,p,u) = ^S(£) 


f 4 

'1 . 

1 

( V \ 

\ - 

- Lm + 0(e) 

+ - 

41n ( 7> ) + 3 

1*7 

e 

e 

V Q ) 


( v\ 9 27r 2 1 

—4I>m !n(— J — 3Lm + 2 g— J . (82) 




O-sCf 

47T 


X 



£ 


e 


+ 


4 Lm + O(e) 
4 L m > ■ 


4 

e 


4 L m 



(83) 


We see that g^ q and are free of large logarithms for 
v = v g ~ Q and v = vs ~ v$ ~ Q(l — x), respectively. 
For later reference we also give the resulting MS-type 
counterterms (subtracting the 1 /e and 1 /77 divergences) 


y<t> 




47T 


1 

+- 

e 


77 [e 


4 In ( 4 


- Lm + 0(e) 


v(i) _ a » C F 

U * Z s ~ — 


Q 


m 


1 


— Lm + 0(e) 


--In 

e 


m 


-i +- 


(84) 


(85) 


B. Two-loop results for the PDF soft and collinear 
functions 


We calculate the secondary massive quark corrections 
at 0(o^CfTf) for the PDF soft and collinear functions in 
the (raj + 1 ) flavor scheme using the (partonic) results for 
the real and virtual radiation of a massive gluon at 0(a s ) 
in Eqs. d and d- Applying dispersion relations as 
discussed in Sec. IV A of Ref. [f| yields 


4/gTF bar<i ’ 2) TO ’ <9- A ’ v) (86) 

= Q> Ab v ) Im [ n (m 2 , M 2 )] 

- ^n( 7 ra 2 , 0 ) - ^ Tf 3 q/g rC l) (^, A, Q, ft, v ), 


S^ +1 ’ h ^ 2 \£,m,A^,F) (87) 

= 5 (bare ’ 1} (^ M, /x, v) Im [n(?7i 2 , M 2 )] 

- (n(m 2 ,0) - a \^ TF ^ 5 (bare ’ x) (£, A, /x, v ), 


with the imaginary part of the vacuum polarization func¬ 
tion n(m 2 ,p 2 ) and its value at zero momentum given by 


Im[ll(m 2 ,p 2 )] = 9(p 

2 3-2d 7 r (3-d)/2 

V 


4 ^) 


d — 


4m 2 )g 2 T F ^(p 2 ) 


4 m* 


) 0 - 


4m 


(d—4)/2 (gg) 

2 \ (d—3)/2 


II(m 2 , 0 ) 




(89) 


For the scheme change contributions in the respective 
second terms of Eqs. ( 15 () 1 ) and ( 1571 ) we use a gluon mass 
A <C 771 as an infrared regulator which allows us to fac¬ 
torize also these corrections with respect to rapidity and 
to use the results from Sec. IV Al The total bare results at 
0(ci 2 CfTf ) for the partonic collinear and soft functions 
then read 


13 Here the e-dependence in the expression proportional to 1/r, e/l] in the ^anomalous dimension. However, for convenience we 

should be in principle kept unexpanded to avoid terms going like show Qnly the terms up to 
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~(n7+l,bare,2) 

U 9 q/q ,T F 


(ai ni+1) ) 2 C F T F _ f 1 


167T 2 


vS. 


(n; + l,bare,2) 
<]>,T f 


1 », ( v 
e 2 [3 n VQ 

80 224 \ 

9 m + 77 J 

yi ni+1) ) 2 C F T F fi 


1 

+ - 
e 


3e 2 9e 3 Lm A+ 2, Lm 


40 




1 47T 2 


80 T 224 X 

T im + 77 + ° (e) . 

_ 16 , f ,8 2 

^ ~r 


\, fv\ 

2 20 

73 

20t r 2 8 > I 

) 111 \o) 

— 4:L m L\ + 2 L m + —L m + 

18 + 

27 3 C3 J 


167T 2 


V 


m 


8 40 16 r r 8 2 80 r 224 

- n 4“ 4“ L m L A — —Lm — L m — - 4" (41(7 

3e 2 9e 3 3 m 9 27 W 


\m 1 


t 

4- 


40 16 


80 


— ■ 7:- 7 rL m L\ + —L + -—L m + 


3e 2 9e 


9 


224 

77 


(90) 


(91) 


with £ = t/v, L m = In(m 2 //r 2 ) and L\ = ln(A 2 //i 2 ). Therefore, the nonvanishing two-loop counterterm contributions 
in MS-renormalization (subtracting also the 1 /rj- divergences), which are used in the (rq + 1) flavor scheme above the 
quark mass threshold, reacF^l 


v z (ni+i,2) = a { r i+1 \^)a { r i+1) ^)C F T F §{ 


j 9,Tf 


In , 

3 V Q 


167T 2 

V 


r 1 r 8 40 16 r 8 o 80 

s — —f — — — —-(- —L T — L ri 
(Tj l 3e 2 9e 3 3 m 9 


224 

77 


0(e) 


-—In 
9 V Q 


1 47T 2 
3 9~ 


(92) 


vZ 


(n,+1,2) _ a 7 ,+1) (/r, 7 «i ra,+1) (M)gF^ / 1 


S,7> 


0 ® 


167T 2 

' 8 40 

37 ~~ 97 


-m 

V 


8 40 16 r r 8 r2 80 r 224 ^,4 

" 37 + 9 ^ + Y LmLA ' 3 Lm ~ y Lm _ 77 + ° {e) 


(93) 


The sum of the individual counterterm contributions 
gives the complete PDF counterterm at 0(a 2 ) with re¬ 
spect to one flavor, 


z (n,+1 ,2) f £_ 

*’ Tf \Q 


1 



y(ni-\- 1 , 2 ) 


y{ n l~\~ 15 2) 

^ S,T f 


) 


(a 7 i+ 1 ) ) 2 C f T F 

167T 2 


|<5(1 - z) 


1 41 



+ 


71-7 

1 - 2 


3e 2 



(94) 


This yields, using also the corresponding contributions 
at 0(a s ) in Eqs. (j84j) and (l85l) the correct 0(a 2 C F T F ) 


contribution to the /r-anomalous dimension for the PDF, 


7&f I,21 (l-z,M) = 

'4 167T 2 


(af l+ ] ) C f T f 
167T 2 

71-7) • 


2 r 


( 2 ) 

Tf 


fl(i-7' 


1-z 
(95) 


-i + 


with = —80/9 being the 0(a 2 C F T F ) coefficient of 

the cusp anomalous dimension Tcusp- 

Below the mass threshold in the (rq) flavor scheme 
the OS subtraction prescription is employed for both the 
strong coupling and the massive quark contribution to 
the collinear and soft PDF functions. The OS prescrip¬ 
tion implies that the secondary massive quark corrections 
decouple in the limit m —> 00 . Since the bare result given 
in Eq. (1901) agrees with its large mass limit it can be easily 
seen that the 0(a 2 C F T F ) massive quark corrections are 
subtracted away entirely by the OS counterterm. For the 
PDF threshold corrections at 0{a 2 C F T F ) we therefore 
obtain 
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We indicate explicitly that only the strong coupling related to 
the interactions of the gluon to the primary quarks is affected by 
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M^(£, m, Q, fj m , v) = J &£' g { ™ l + 1] (£-£', m,Q,k,ii m ,v) 


O(al) 


^(nj+l,bare,2) ^(n^ + 1,2) -f .F 


9q/q,Tp 


9,Tf 


37T 




/ 4bare,l) 7 ( ni +l,l)\ 
\ 9 q/q ~ X 9 J 


A4g 2 ) (£,m, = f d £>S^\£-£>,m,A,, m ,v)(s^) (£', A, /j m , v) 


o{ a f, 


(rai+l)™ 

_ o(^i + libare,2) _ ^(ni + 1,2) _ ^^ j { Afbare,!) r^Gi + l,!) 


= 5, 


Tp 


z : 


s,t f 


37T 


A, 




>)■ 


(96) 


(97) 


Note that the difference of the scheme for a s in 
g^JJq, and g^ +1 \ 5C-+ 1 ) affects the terms at 

0(oi 2 s CfTf ) and leads to the third term in the last equal¬ 
ity of Eqs. (1961) and (ED, respectively. Using the two- 
loop results in Eqs. (liffll) (l93l) and the one-loop results in 
Eqs. (1821) (l85l) with a gluon mass A as an infrared regu¬ 
lator we obtain 


Vet, M {2) (£,m, = 


(a { r i+1) ) 2 C F T F 


8 2 80 r 
+ -rrL r 


3 m 9 

2 8tt 2 

3 + ~9~ 
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167T 2 
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224 

~TJ 


(99) 


where the dependence on the IR regulator A has dropped 

( 2 ) ( 2 ) 

out. Upon summing up M g and Mg we obtain the 
total PDF threshold correction already given in Eq. ED- 
The RGE for the ^-evolution of the threshold correc¬ 
tions reads ( i = g , S) 


v-j- Mi(e,m,Q,/jL,u) = ^Mi Mi(£, m, Q, fi, v ). (100) 
d v 

The jx-anomalous dimensions 7 a^ s and iMs can be di¬ 
rectly read off from Eqs. (l98l) and (l99l) or equivalently 
from the ratio of the counterterms in the (n{) and (rq + 1 ) 
scheme for g and S in analogy to the last equality in 


the rapidity regularization procedure and adopts a ^-dependence. 
The interactions due to gluon splitting within a single sector do 
not contain any rapidity divergences and therefore do not need 
additional regularization. We note that the renormalized strong 
coupling depends on the scale v only due to the dimensional 
extension of the k ~-integration and satisfies da s (/x, v)/&\nv = 
— r]Oi s (ii, u) to all orders. 


Eq. ED- which gives 


iMg = -7 Ms 

01 2 s CfTf 

167t 2 




+ 0(a 3 s ). 

( 101 ) 


The solution of the rapidity RGE in Eq. (11001) is a simple 
exponentiation of the rapidity logarithm, i.e. setting v = 
v g in A i g and v = = us in M.$ we get 


- z, to , fj, m ) = Q J d£Ms{£,m,iJ, m ,vs) 
x M g (Q(l — z) — £, Q, to, /i m , v g ) (^“^ • ( 102 ) 

In order to allow for an arbitrary evolution path in /x- 
TOspace one can generalize this expression to resum the 
logarithms In (m 2 /g) in the i/-anomalous dimensio n by 
integrating the latter in /x as discussed in Refs, fill I23|. 
The variations of the scales us and u g may be used as 
an additional input for the perturbative uncertainty es¬ 
timate. Finally, we remark that by setting v g = Q and 
vs = u,f, one can obtain the compact all-order expression 


AdJ(l - z,m,nm) 
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IMg MfaS , 

(103) 


where 7 x s is the ^-anomalous dimension with the 2 - 
loop contribution given in Eq. (11011) . denotes the 

coefficient of the ^-distributions in the PDF threshold 
correction, i.e. 


■M.Qfi — 1 + 


cx 2 CfTf 
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(104) 


The noninteger plus-distribution in Eq. (11031) is defined 
as the analytic continuation of 0(1 — z)/( 1 — z) 1 ~ 1M n , see 
the appendix of Ref. [ 6 ] for details. Expanding Eq. (11031) 
in a s allows one to easily to read off the distributive 
structure of A4j at any order in the strong coupling in 
terms of the anomalous dimension 7 _m 9 ■ 
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C. Threshold corrections for N 3 LL analysis 


For a complete analysis at N 3 LL we need the terms at 
0 (o^ln 2 (l — x)) and 0 (a 3 ln(l — x)) both counting as 
0(a 2 ) for a s ln(l — x) ~ 1. The former can be easily ob¬ 
tained from the exponentiation property of the rapidity 
logarithm. The latter can be read off from the nonsin¬ 
glet PDF threshold correction in the OPE region that has 


been recently computed up to 0(a: 3 ) in Ref. [25j]. The 
corresponding expanded result for x —> 1 (Eq. (5.60) in 
Ref. [25|) fully agrees with our computation for the /z m - 
dependent terms at 0 (a 3 ln(l — x)) (which are obtained 
from the ratio of the evolution factors in the ( iy ) and 
(ni + 1 ) flavor schemes), but in addition allows us to ex¬ 
tract the relevant /i m -independent term. 

The complete result for the PDF threshold correction 
in the logarithmic counting a s ln(l — x) ~ 1 at N 3 LL 
reads, (£ = £/us, £ ~ vs ~ ~ <9(1 - x), v g ~ Q) 
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■ 77 -Mil --,m, Q,n m ,vg,vs ) = S(£) + 
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J 0(a 2 ) 


where the second term in the first line counts formally as 0(a s ) and is therefore already relevant at N 2 LL. The 
(universal) functions related to the rapidity logarithms read (L m = ln(m 2 //x^)) 


A41 ln( m ’ — CfTf -F^rn 4-+ 


l 0,ln 


{m,Hm) =C f Tf <{ L 


80 224) 

— . _i- K 


9 


27 J 


(106) 


l 0,ln 


176 64 128 ' 

-vr CA + V7 TFni+ vF TF 


T Lm 

+ C A 


/1240 160tt 2 224 

\~81 27 + ~~3 

/ 35452 1648tt 2 

y - 729 243 


+ K 


4 —— C3 ) Ca 4- ( — — 64^3 ) Cf 




-24 C f +™T f 


2176 1984 

■ Tpni H——— Tp 


81 


81 


-6OC3 


176t r 4 32 


where 


/24064 512 / 

V 729 ~~ 

( M< £}n (m,Ai m )) 


+ Tpni ^ 72Q 27 ( 3 J + T Fy 72g 4- ^ Cs 


135 3 

12064 896 


-7T-® 4 ) 4 - Cp ( — 


2834 1208 . 16tt 4 64 


27 


C3 ^ 


15 


T* 4 


32 


640 


= C Z pT^—Lt + — L 


1664 


27 


27 


■ L rn + 


17920 


243 


25088) 
729 J 


B 4 = - In 4 (2) - — In 2 (2) - 

3 w 3 w 180 


+ 16 


u,Q. 


The PDF specific function at 0(a 2 ), which is not multiplied by a rapidity log, reads 
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For completeness we display also the threshold corrections for the hard and jet functions. As discussed in Ref. [g] 
the characteristic rapidity scales for the mass shell fluctuations in are v 4 ~ Q and ~ m (with the 

symmetric tj- regulator), while in A4y they are ~ Q and jyf ~ m 2 / (<9(1 — a:)). To account for correlations between 
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these scales we set = v{ = z / 9 and 1/3 = v^/vg. With these choices the results at N 3 LL have the form 
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As stated at the end of Sec. m the consistency rela¬ 
tion (l62l) implies that the coefficients of the rapidity log¬ 
arithms, i.e. the ^-anomalous dimensions, are the same 
for all threshold corrections, i.e. 

■ M H ) ,ln( rn ^m) = Hm) = fX m ) , 

(113) 

■Wffint" 1 -/ 1 ™) = M Jln K/ 1 ™) = M ^\n 

(114) 

Mglnifw^m) = M { ^ ln2 {m, Hm) = Hm) , 

(115) 


which has been already used implicitly in Ref. [fj. Fi¬ 
nally, the remaining function-specific 0(a 2 ) corrections 
read 
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Since hard and jet functions and PDFs are building 
blocks of factorization theorems for many different pro¬ 
cesses at hadron-hadron collisions, the results for the 
threshold corrections can be directly applied there as 
well. 


VI. CONCLUSIONS 

In this work we have discussed how to set up a VFNS 
for multiscale processes at hadronic collisions, where we 
have taken inclusive DIS in the endpoint region x —> 1 
as a specific example. In this limit massive quarks do 
not (predominantly) participate directly in the hard in¬ 
teraction with the virtual photon and therefore mainly 
arise as secondary radiation giving corrections starting 
at 0(a 2 ) in the fixed-order expansion. Starting from the 
massless factorization theorem we have shown how to 
systematically incorporate the secondary massive quark 
effects by using two kinds of renormalization conditions 
for the massive quark corrections in the gauge invariant 
components. The use of the MS renormalization pre¬ 
scription in the small mass region and of the on-shell 
(low momentum subtraction) renormalization prescrip¬ 
tion in the large mass region imply automatically that 
all large logarithms are resummed and that the respective 
correct limiting behavior is achieved in the massless and 
the decoupling regions. The difference between these two 
schemes manifests itself in additional threshold matching 
corrections at the mass scale. We have discussed some 
universal features of these threshold corrections, which 
exhibit intrinsic relations among each other due to the 
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consistency of RG running. Here we have also computed 
explicitly the PDF threshold correction for x —> 1 at 
0(a 2 ) and showed how to resum a large remaining log¬ 
arithm therein that is related to the separation of mass 
shell fluctuations along rapidity and displayed final ex¬ 
pressions for a N 3 LL analysis. From a practical point 
of view our VFNS in the endpoint region of DIS can be 
combined with a VFNS in the OPE region 1 — x ~ 0(1) 
by adding the known associated nonsingular corrections 
related to the difference between the full perturbative 
QCD result and the fixed-order expressions for the com¬ 
ponents of the SCET factorization theorem for x —> 1. 
This may have an effect also for moderate values of x due 
to dynamical threshold enhancement (see e.g. Ref. M), 
an effect which reinforces perturbative corrections close 
to the partonic threshold due to the steep fall-off of the 
PDFs for momentum fractions close to one. 

While we have concentrated in this work on DIS, the 
concept of how to theoretically treat the effects of sec¬ 
ondary massive quarks within factorization is applicable 
for more general processes including hadron-hadron col¬ 
lisions. In particular, the massive quark threshold cor¬ 
rections relevant for resummation of logarithms at N 3 LL 
order determined for the massive components in the fac¬ 
torization theorem (hard function, jet function, PDF) are 
universal and can be employed in factorization theorems 
for other processes where these components appear. 


where the index q includes both quarks and antiquarks. 
The factorization theorem for F 2 (x, Q) is analogous. 
Note, however, that the Callan-Gross relation F 2 = 2xF\ 
does not hold in the OPE region beyond tree level. In 
the Breit frame the PDFs fare forward matrix el- 

J o/p 

ements of SCET operators decomposed out of collinear 
fields [13, HI] ■ For definiteness we set the final renormal¬ 
ization scale to be p, = /j,h- In the massive quark case 
one obtains 


d£ fdC 


Fl(x,Q,m)= f Y J c J ci 
i=q,Q j,k=q,g J 


(A2) 


x H, 


(«i) 




for fj, m > hh and 



x Mf,u (^—,m,Urn^j Ufj™ — fh/piZ ■ 


ACKNOWLEDGMENTS 

We would like to thank Ilaria Jemos for collaboration 
at an early stage of this work. The work of P.P. was sup¬ 
ported by the German Science Foundation (DFG) under 
the Collaborative Research Center (SFB) 676 Particles, 
Strings and the Early Universe. D.S. is supported by the 
FWF Austrian Science Fund under the Doctoral Program 
No. W1252-N27 Particles and Interactions. 


Appendix A: Secondary massive quark corrections 
in the OPE region 


for fx m < fin. The secondary massive quark corrections 
to the hard functions and F[^ l+l ' > at 0(a 2 CpTp) 

can be written as 

H^\z, Q,m,n)= H™(z, Q, m) + F^ 2) (z, Q, m ), 

(A4) 

H <” i+1) (z, Q,m,n)= i^ !+1) (z, Q,/x) + P[ n ^ 2) (z, Q, m ). 

(A5) 

The full QCD result at 0(a 2 s CpT F ) is both IR- and UV- 
finite and can be decomposed into a purely virtual cor¬ 
rection and a real radiation correction with the kinematic 
threshold z = 1/(1 + 4m 2 ), 


We display explicit results for the perturbative correc¬ 
tions due to secondary massive quarks at 0(a 2 ) in the 
classical OPE region, where 1 — x ~ 0(1). Here the fa¬ 
miliar factorization theorem for n / massless quarks reads 


*iM) = £| £ 


i=q 


3=9,9 


rl di (n/ ) 
„ £ ij 




(Al) 


F^ 2 \z,Q,m) = 2F^‘’ 2 \Q,m)S(l - z) 

+ 9(z)9( 1 - z - 4m 2 z) (z, Q, m ), (A6) 

where rh = m/Q. The partonic QCD current form factor 
Fm l ' 2 \Q,m) is given in Eq. (l30l) . The real radiation 
function F^ l g 2 \z,Q,m) was first computed in Ref. 0 
(and also checked by us) and reads 
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Here we have used the abbreviations 


r z = y/l — 4:iri 2 z , w z = i/l- 
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We remark that the QCD corrections decouple in 
the heavy quark limit using the ni scheme for a s , 
i.e. F < ^ 2 \z,Q,rn) —> 0 for rh —> 00 . In the small mass 
limit rh —> 0 , on the other hand, we obtain 
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The function represents the quark mass correc¬ 

tion to the massless quark result in the (ni + 1 ) flavor 
scheme and 43 ! 
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The PDF threshold correction Mf qq is given by 
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Here the superscript ni + 1 indicates only that the (nj + l) scheme where the scheme for a s does not need to be specified at 


for Oi s is used in the expressions m and ED . 


this order. 
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FIG. 7. Secondary massive quark contributions at 0(o. 2 CfTf) for rh — m/Q = 0.1 in full QCD (blue, solid) together 
with the singular result for z —> 1 at fixed order (red, dashed) and the nonsingular terms (green, dotted), all normalized by 
a 2 CfTf /(47t) 2 . The left panel shows the purely partonic result with initial state quarks, while in the right panel we convoluted 
the partonic form factors with the function (1 — x) A representing the steep decrease of the PDFs for large values of x. 


Appendix B: Expansion for x —¥ 1 

The massive quark corrections to the factorization 
theorem discussed in Sec. um represent the singular 
0(oqCfTf ) secondary massive quark corrections to the 
structure function F\(x,Q,m) in the fixed-order expan¬ 
sion in full QCD. Besides the virtual contributions in 
QCD, which are fully contained in the SCET description, 
the singular perturbative fixed-order corrections also con¬ 
sist of the collinear real radiation contributions which 
arise for 1—x ~ m 2 /Q 2 <C 1. Setting p, = /ih = ftj = /x m 
in Eq. (l25l) we obtain 

F 1 (x,Q,m)\ FO = ^2^- f&Zcttyp - x, m ) 

i=q 

x Q 2 H ( ' ni \Q,m,ii) J (ni) (Q 2 (l - . (Bl) 

'-V.-' 

where the massive quark contributions to the fixed-order 
hard function at 0(cx 2 CfTf ) read 

=2F^ 2 \Q,m)6(l - z) 

+ Q 2 J^{Q 2 (l-z),m). (B2) 

with Fjn l ’ 2 \Q,m) and ^^(s,m) given in Eqs. (l30l) 

and C!?l) with a s = or™ 1 '(ft). We can obtain this result 
also from the corresponding full QCD fixed-order result 
in Appendix[A] Since the virtual contributions multiplied 
by the <5(1 — z ) distribution in the OPE and endpoint 
regions agree, we only have to consider the expansion of 


the real radiation term F^ n // 2 J (z , Q , m) in Eq. K71) for 
1 — z ~ to 2 < 1 (with m = m/Q), which yields indeed 
the correct term, 

0(1-z- 4 rh 2 z) (z,Q,m) (B3) 

^ Q 2 J^(Q 2 (1 - z),m) + 0(( 1 - z)°, m°). 

In Fig. [7] we investigate how well the expansions work for 
the specific scale ratio to = 0.1. The left panel shows 
the partonic result for the full QCD corrections (blue, 
solid), i.e. F[ ni ' 2 \z 1 Q,m) in Eq. (IA6I) . the singular re¬ 
sult for the endpoint region at fixed order (red, dashed), 
i.e. H^/ I //, 2 \ l rn (z, Q , to) in Eq. (1B2I) . and the difference de¬ 
scribing the nonsingular corrections (green, dotted). We 
see that for 1 — z < 4 m 2 /Q 2 the endpoint corrections 
encode the dominant behavior, but fail to give a good 
description of the full QCD form factor below a certain 
value, here z < 0.5. The right panel displays the ab¬ 
solute value of the convolution between these partonic 
functions and a common function f(x) = (1 — x ) 4 acting 
as a dummy PDF at the endpoint which falls off steeply 
for x —> 1. The convoluted results for the full QCD and 
singular terms are negative for x > 0.05. We see that at 
this level the agreement between the singular and full re¬ 
sults is much better due to dynamical threshold enhance¬ 
ment (see e.g. Ref.[2f|), up to values significantly lower 
than x = 0.5. This may have the consequence that end¬ 
point region effects can have an impact even at smaller 
values of x probed at hadron-hadron colliders. A recent 
analysis on this issue was carried out in Ref. [29}. They 
found that the effect can be sizable and may require the 
use of resummed PDFs for resummed calculations. 
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